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The phase diagram of adsorbate-induced row-type-alignments, such as missing-row reconstructions 
induced by adsorbate-atoms on the FCC (110) surface, is calculated by the Blume-Emmery-Griffiths 
(BEG) model. In the model, we introduce adatom-adatom and dipole-dipole interactions between 
nearest-neighbor (NN) and next-nearest-neighbor (NNN) rows. The calculation of the temperature 
versus adatom chemical potential phase diagram is performed using mean-field approximation. It is 
indicated that when NN and NNN interactions are competitive, there appear either dipole or cover- 
age modulated (incommensurate) phases at high temperatures for wide regime of the interactions. 



I. INTRODUCTION 

It is well known that row-type-alignments are induced 
in many cases on transition metal surfaces by adsor- 
bates. A typical example is the (1 x n) missing- row 
reconstruction (MRR) of FCC(llO) metal surface coh- 
ered by adatoms- such as H/Cu(110)H, O/Rh(110)tH2, 
and O/Pd(110)Dtl. In these systems, the rows of the 
the metal atoms in [110] directions are accompanied by 
zigzag chains of adatoms sitting in the threefold FCC 
hollow sites on the [111] microfacets of the surface and 
each nth row is missing. It should be noted the arrange- 
ment of the missing row at low coverage c of adatoms 
suggests that the interactions between the rows of the 
metal atoms are of repulsive. _ 

On the theoretical side, P. J. Kundrotas et alU re- 
cently showed that the 2 dimensional Blume-Emmery- 
Griffiths(BEG)-like modelcl is a useful model to analyze 
these ordering phenomena of row- type-alignments includ- 
ing the adatom-metal dipoles. Their excellent calculation 
of the finite temperature phase diagram for O/Rh(110) 
is, however, restricted to the case where the nearest- 
neighbor (NN) interactions along [001] direction is much 
stronger than the next-nearest-neighbor (NNN) interac- 
tions. 

The purpose of this paper is to extend the analysis, 
including the ground state analysis, to the case of the 
relatively strong NNN interactions. In the case of axial- 
next-nearest-neighbor-Ising (ANNNI) modelO, when NN 
and NNN interactions are competitive, an introduction 
of the relatively strong NNN interactions usually give rise 
to change of the ground state structures and also induces 
entropic modulated phases at high temperatures. From 
this fact, we may expect that the BEG model will also 
exhibit a rich variety of phases due to the frustration of 
the interactions. 

In this paper, after establishing our model, we deter- 
mine the ground state phase diagram by energy compar- 
ison method. This ground state analysis shows the new 



structure corresponding to the (uudd) structure of the 
ANNNI model appears under the certain conditions for 
NNN interactions. We then derive the mean-field equa- 
tions to investigate qualitative feature of the finite tem- 
perature phase diagrams. Solving the linearized mean- 
field equations analytically, we obtain the condition for 
the interactions to realize the stable modulated phases. 
Numerical analysis of the mean-field equations is per- 
formed to calculate the finite temperature phase dia- 
grams for several interaction regimes, including the re- 
gion satisfying the modulated phase conditions. These 
analyses show that there appear either dipole or cover- 
age modulated phases at high temperatures for certain 
regime of the interactions. 

Organization of this paper is as follows: In the next 
section (section II ), after setting up the model for the 
MRR ordering, we construct the ground state phase dia- 
gram. Section III deals with derivation of the mean-field 
equations and solution corresponding to the modulated 
phase. Section IV is devoted to the numerical analysis. 
Finally, summary and concluding remarks are given in 
section V . 



II. MODEL AND GROUND STATE ANALYSIS 

It is plausible to expect that the clusterlike bonding be- 
tween metal atom and adatom produces a dipole moment 
by transfer of electron. In case of zigzag adsorption, there 
exists the asymmetric bonding between metal atom and 
adatom. The asymmetry can be theoretically treated by 
the dipole moments with opposite directions, say "up" 
or "down" dipoles. Therefore, a site of the surface where 
metal atom sits in is in one of three states, missing of 
metal atom (vacant), occupied by metal atom with "up" 
or "down" dipoles. 

In order to analyze the ordering phenomena of row- 
type-alignments including the adatom-metal dipoles, we 
introduce pseudo "spin" variable Cj at i site to represent 
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the three states of the site i, vacant (cr^ = 0), occupied 
by metal atom with "up" or "down" dipoles (cr, = ±1). 
Since the occupation number of metal atom (accompa- 
nied by adatom) at i site is given by of, the Hamiltonian 
of the system is represented by the BEG model as 



H = ^JyWj + - 2 / 3 )K 2 - 2 / 3 ) 



2/3) 



(1) 



where Jij (Kij) 



is interaction between ij pair of the 
dipoles (metal atoms), /i is the chemical potential. <rf — 
1(0) means a metal atom occupies the surface site i (or 
not) . Here we use (of — 2/3) instead of of for convenience 
of the later calculation. Following P. J. Kundrotas et aM, 
we assume that 



(Jij 7 Ki : 



(Jo,Kq) ij = NN pair in [110] direction 
(Ji,Ki) ij — NN pair in [001] direction 
(J 2 ,K 2 ) ij = NNN pair in [001] direction 
(0, 0) otherwise 

(2) 



and farther neighbor interactions are negligible. Then 
the Hamiltonian is rewritten as 

U = ^ j J Q a em oe m+i + Ko(<rj m - 2/3)(cr 2 m+1 - 2/3) 
+ J\<?lm&l+\m + Ki( a £ m ~ 2/3)(cr 2 +lm - 2/3) 



+ M (<rf ro - 2/3) j 



2/3)K 2 +2m -2/3) 



(3) 



where I (m) refers site-number along [001] ( [110]), 
namely I- (m-), direction. 

In the MRR-type-alignments, metal atoms form a 
chain-like structure and the adatoms are ordered in 
zigzag arrangement which corresponds to " antiferromag- 
netic" ordering of the <jg rn along m-direction. From these 
experimental facts, the interactions in Eq.(||) should sat- 
isfy the 



J > , K Q < , Ki > , K 2 > 



(4) 



It should be noted that the Hamiltonian (|J) can be trans- 
formed into "ferromagnetic" (Jo < 0) one by the trans- 
formation a im — > (-l) m (T£ m . 

Now we determine the ground state phase diagram by 
comparing the energies of several "spin" arrangements. 
From the condition J > and Kq < , the "spin" 
arrangement in the ground state is expected to be repre- 
sented by m-independent "spin" &i as 



&tm = (— l) m <5^- 



(5) 



To represent the "spin" arrangement, we introduce the 
notation which indicates the repeating unit of og as 
(<7iCT"2 • ■ • (J n ); for example (udO) represents the "spin" 
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FIG. 1. Ground state phase diagrams for Ko < 0, 
Jo > and Ki > |Jj|(i = 1,2) in k 2 -\i 2 plane, where 
k 2 ={K 2 - | J 2 |)/(A"i - |Ji|) and ^={-^+±{K + K x + K 2 ) 
+ \Jo\ - Kq)/{Ki - \ Ji\): (a:upper) for Ji > 0, J 2 > 0; 
(b:lower) for Ji < 0, J2 < 0. Boundary lines in the diagrams 
are described as: fi 2 = for (0) and (ttOO) ; fJ, 2 = 3k 2 for (uOO) 
and {uOdO) in (a), for (u00) and (u0) in (b); fj, 2 = 3/2 for (u00) 
and (ttdOO) in (a), for (ttOO) and (uuO) in (b); fj, 2 — 2 — k 2 
for (uOdO) and {udO) in (a), for (uO) and (uuO) in (b); 
k 2 = fJ. 2 /2 — 1/4 for (tidOO) and {udO} in (a), for (uuOO) and 
(uuO) in (b); = (fJ, 2 — Aj)/2 for (WO) and (ud) or (uudd) 

in (a) where A, = Min(%.hl(K x - \ J x |), 3Ji/(iTi - |Ji|)); 
&2 = (a*2 ~ 2)/2 for (umO) and (u) in (b). 

arrangement with repeating unit ((Tito's) = (+1 —10). 
The target " spin" arrangements of the energy comparison 
are 

Case A: Ji > , J 2 > 
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(0) , (uOO) , (uOdO) , (udOO) , (udO) , (ud) , (uudd) 
Case B: Ji < , J 2 < 

(0) , (uOO) , (uO) , (uuOO) , (uuO) , (u) 
It is noticed that the results for the case C : J\ < , J 2 > 

(case D :Ji > , J 2 < ) can be easily obtained 
from the results for case A (case B) by the transformation 
at — » (-1) £ <7£. 

Figure 1 shows the ground state phase diagram in u 2 - 
k 2 plane obtained under the condition Ki > \Ji\(i = 
1,2), where fi 2 and k 2 are given as 

k 2 = (K 2 -|J 2 |)/(^i-|Ji|), (6) 

M2 = (-^+f(^0+^l+^ 2 ) + 

|Jo|-^ )/(^i-|Ji|). (7) 

For large /i 2 regime in Fig. 1(a), either (ud) or (uudd) is 
realized depending upon 2J 2 < J\ or 2J 2 > Ji, respec- 
tively. 

For large K 2 , the energetically preferable structure 
may be the one containing no K 2 bonding, such as (uuOO) 
or (udOO). Actually, as is shown in Fig. 1, there appears 
either (uuOO) or (udOO) depending sign of J\ for the rela- 
tively large K 2 . Therefore we might observe the different 
sequences of the MRR with increasing the chemical po- 
tential if we prepare the surface with relatively large K 2 . 
It is remarked here that there exist many structures de- 
generate in energy in the region of the phase containing 
two adjacent "0". For example, structures consisting of 
any combination of (uOO) ((udOO)) and (dOO) ((duOO)) are 
degenerate in energy with the structure (uOO) ((udOO)). 

III. MEAN-FIELD ANALYSIS AND 
MODULATED PHASES 

In this section we derive the mean-field equations to 
construct phase diagrams at finite temperatures. The 
mean-field energy e per m-direction line is obtained as 

e= ^{-|J |^ + ^( y | + 2 ) 

+ J&txt+i + ^-(ym+i + 2) 

+ J 2 X£x e+2 + ^-(y e y e+2 + 2) - (8) 

where xe and ye are m independent thermal average of 
(— l) m ai m and (2 — 3erf TO ), respectively. Since af m is 
the occupation number of metal atom (accompanied by 
adatom) at Im site, as mentioned in previous section, 
ct = (2 — yi)/3 gives the m independent thermal average 
of atomic concentration of metal or adatom (coverage) at 

1 m site. 

The mean-field entropy s per m-direction line can be 
expressed in a unit of ks = 1 by xt and yi as 

s = -££ik(a){]ogik(a) - 1} . (9) 

t cr 



where Rt(a) is the m independent probability of finding 
the £m site at state a given as 

Rt(a) = ^(1 + 3axt/2 + (2 - 3<7 2 )^/2). (10) 

Optimization of the free energy / = e — Ts per line 
with respect to xg and ye leads the so called mean-field 
equations as 

- 2| J \x e + Ji(xt-i + xt+\) + J 2 (xe- 2 +xe +2 ) 

+ T^|lo gj R,(a)-0, (11) 

CT 

2-^o K lf , K 2 , 

—Q-vt + -g-w-i + Vt+V + "g-W-2 + 2^+2 + 2) 

- f + T^^^log W =0. (12) 

where T is the temperature. 

When a = 0, we easily obtain a stable "paramagnetic" 
solution xe = and ye = with coverage c = 2/3 at 
high temperatures. To determine the transition temper- 
ature from the "paramagnetic" phase and the phase just 
below the transition temperature, we linearize the mean- 
field equations. Substitution of xe — X q exp(iq£) and 
ye = Y q cxp(zg£) in the linearize equations leads us to 
the transition temperature at a = as maximum of 

T x {q) = |(|Jo| - JiCos( g ) - J 2 cos(2 g )) , (13) 

T y (q) = 1(-K - K lC os(q) - K 2 cos(2q)) . (14) 

When NNN interactions J 2 | , \K 2 \ are very small, ei- 
ther q = or q — tt gives the maximum and there ap- 
pears no modulated phases. On the other hand, if NNN 
interactions satisfy the condition 

J 2 /Ji| > 1/4(J 2 > 0) or \K 2 jK x I > l/4(^ 2 > 0) 

(15) 

the system just below the transition temperature ex- 
hibits the modulated phase characterize by wave num- 
ber cos(Q) = — J1/4J2 (dipole modulation) or cos(Q) = 
—^1/4^2 (coverage modulation) . 

IV. NUMERICAL RESULTS 

Since it is difficult to solve the mean-field equations an- 
alytically, we have to solve the equations numerically to 
obtain the phase diagrams. In the present paper, we con- 
fine ourselves to the case of " ferromagnetic" interactions 
between rows, i.e J\ < and J 2 < 0. To demonstrate 
the change of ground state and appearance of the mod- 
ulated phase, we deal with the three cases of interaction 
regime, 

a) J = 1, K = -3, Ji = -1, K x = 3, J 2 = -0.05, 
and K 2 = 0.5 , 
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-0.05, 



Jo = -0.03, 



b) Jo = 1, K 
and K 2 = 1.55 , 

c) J = 1, K Q = -3, Ji = 
and X 2 = 0.88 . 

The interaction constants in the case a) are the same 
as Rcffl; those in the case b) correspond to the regime 
k 2 = 0.75(> .5) where the (uuOO) appears as ground 
state (see Fig. 1(b)); those in the case c) satisfy the condi- 
tion K 2 /K t {= 0.29333) > 1/4(K 2 > 0), the modulated 
phase is expected to exist. 

To construct the phase diagrams for the three cases 
of interaction regime at finite temperatures, we solve 
the mean-field equations within repeating unit 12, i.e 
(%i2+e,yi2+e) = {%t,Di) {i= 1,2, ■•■,12), for various 
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FIG. 2. T versus fi phase diagrams: (a:left upper) for 
Jo = 1, K = -3, Ji = -1, Ki = 3, J 2 = -0.05, and 
K 2 = 0.5 ; (bdeft lower) for J = 1, K = -3, Ji = -I, 
A'i = 3, J 2 = -0.05, and K 2 = 1.55 ; (crright) for J = 1, 
Ko = -3, Ji = -0.5, K t = 3, J 2 = -0.03, and K 2 = 0.88 . 
The dotted(full) lines represent the 2nd(lst) order phase tran- 
sition lines. 

values of temperatures and chemical potential /12 in 
Fig. 1(b) along the line fc 2 = 0.225 ( case a)), the line 
k 2 = 0.75 (case b)) and the line k 2 = 0.34 (case c)). 

The calculated T versus \i phase diagrams for the cases 
a), b), and c) are shown in Fig. 2(a), (b), and (c), respec- 
tively. In the figures, u/d denotes the "paramagnetic" 
row with x = and t; / 0. 

It is natural to expect that at low temperatures, 
there appear successively the phases corresponding to the 
ground state structures with increasing chemical poten- 
tial; at high temperatures the rows prefer the " paramag- 
netic" state to "magnetic" state to gain free energy by 
entropy term. As seen in Figs. 2(a) to (c), all of our phase 
diagrams show this tendency. 

It is clear that the phase diagram in Fig. 2(a) is the 
topologically same as that in Ref.Q. This confirms qual- 
itative feature of our calculation. Major differences be- 
tween our results and those of RefE are the transition 
temperature of our mean-field calculation being as much 
as twice of that in Ref.Q and existence of the " magnetic" 
phase (uOO) phase at low temperatures. 

As shown in Fig. 2(b), (uuOO) phase corresponding to 
the ground state structure in which two out of four rows 
of metal atoms are missing appears for a certain range 
of chemical potential at low temperatures. P. J. Kundro- 
tas et alU did not take this (uuOO) into account. From 
predictions, the modulated phase accompanied by the 
wave number Q — cos -1 (— Ki/4K 2 ) should appear in 
this case. However, we could not find it because our 
calculation is limited within the periodicity 12 and not 
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enough big to find it. 

Fig. 2(c) shows the phase diagram for the case c) where 
the condition K 2 /Ki > 1/4 (K 2 > 0) is satisfied. There- 
fore, the modulated phase with the wave number Q — 
cos -1 (—K1/4K2) — Qi = 107r/12 should appear in this 
case. Actually, we obtain within repeating unit 12 the 
stable modulated phase shown in Fig. 3 in the vicinity 
of the point (// = , T y (Qi)). In the figure, r repre- 
sents the row numbers along the [001] direction on the 
FCC(llO) surface and y r = 2 — 3c r where c r describes 
the probability that a metal atom row exists at the rth 
row. In comparison with the structure (uuO) at lower 
temperatures in the case (shown in Fig.2(c)), the dis- 
tribution of metal atom rows along the [001] direction 
on the FCC(llO) metal surface is modulated. Thus our 
numerical results are consistent with the predictions. Al- 
though it will be difficult to find out experimentally the 
modulated phase, we show a possibility of realization of 
the modulated phases. 

The details of the present calculation will be published 
near future 

1, , , , , , , , , , , , 




appears under the certain conditions for NNN interac- 
tions. 

For the finite temperature, we obtain the condition for 
the interactions to realize the stable modulated phases by 
solving the linearized mean-field equations analytically. 

Numerical analysis of the mean-field equations is per- 
formed to calculate the finite temperature phase dia- 
grams for several interaction regimes, including the re- 
gion satisfying the modulated phase conditions. These 
analyses show that there appear coverage modulated (in- 
commensurate) phases at high temperatures for certain 
regime of the interactions. 
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FIG. 3. r-dependence of y r in the modulated phase. The 
dotted line represents 0.6 cos(Qi(r — 6)) — 0.1 cos(Q2('" — 6)), 
where Qi = lOvr/12, and Q 2 = 9tt/12. 



V. SUMMARY 



The ground state and finite temperature phase dia- 
grams of adsorbate-induced row-type-alignments is cal- 
culated by the Blume-Emmery-Grifhths (BEG) model. 
In the model, we introduce adatom-adatom and dipole- 
dipole interactions between nearest-neighbor (NN) and 
next-nearest-neighbor (NNN) rows. 

The ground state analysis shows the new structure cor- 
responding to the (uudd) structure of the ANNNI model 
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